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Abstract. We investigate Lemaitre-Tolman-Bondi (LTB) models, whose early time evolu¬ 
tion and bang time are homogeneous and the distance - redshift relation and local Hubble 
parameter are inherited from the ACDM model. We show that the obtained LTB models and 
the ACDM model predict different relative local expansion rates and that the Hubble func¬ 
tions of the models diverge increasingly with redshift. The LTB models show tension between 
low redshift baryon acoustic oscillation and supernova observations and including Lyman-a 
forest or cosmic microwave background observations only accentuates the better fit of the 
ACDM model compared to the LTB model. The result indicates that additional degrees of 
freedom are needed to explain the observations, for example by renouncing spherical sym¬ 
metry, homogeneous bang time, negligible effects of pressure, or the early time homogeneity 
assumption. 
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1 Introduction 

Supernova la (SNIa) observations at the end of the 1990’s contain the peculiar feature that 
the universe appears to be accelerating [1, 2]. In spite of years passing by, the phenomenon 
remains without complete understanding; the issue is generally recognised as the dark en¬ 
ergy problem. Attempts to solve this problem include modifying gravity and relaxing the 
cosmological principle, but so far has no single theory explained the accelerating expansion 
together with the other cosmological phenomena better than the ACDM model. 

In the framework of Einstein’s general relativity, a natural way to create apparent accel¬ 
eration is to use large-scale inhomogeneous matter distributions [3], which means giving up 
the cosmological principle. As the principle is not observationally conhrmed [4, 5], renounc¬ 
ing it is merely a philosophical question and as such a widely considered option. However, 
to describe an inhomogeneous matter distribution in general relativity is mathematically 
a very challenging task, e.g., the linear perturbation theory of the Lemaitre-Tolman-Bondi 
(LTB) models [6-8] is not completely understood [9], even though the LTB models are con¬ 
sidered as simple inhomogeneous models. This challenge is to be taken step by step, revealing 
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the possibilities brought by inhomogeneously distributed matter starting from the simplest 
models. 

A vastly investigated LTB model is where the observer is located inside a void. It is 
often assumed that the void is embedded in a homogeneous and isotropic background, and 
that the void was shallow enough to be indistinguishable from its background at early times. 
This simplification allows one to apply homogeneous and isotropic perturbation theory. The 
observed isotropy of the cosmic microwave background (CMB) enforces the observer to be 
located very close to the center of the void [10, 11]. Attempts have been made to find models 
where this property is removed, but it appears that the LTB models do not have enough 
degrees of freedom to allow it [12]. 

LTB models, where the observer is located at the symmetry center, are widely studied. 
The authors of Ref. [13] found SNIa, first CMB peak and baryon acoustic oscillation (BAO) 
observations to fit the LTB model they used. In Ref. [14], their model was able to provide 
good fits to the same kind of data, but encountered difficulties with the local value of the 
Hubble parameter (Hq). In the models used in Ref. [15], the authors found strong tension 
between SNIa and BAO observations, which became even worse when CMB observations 
was included. Additional problems were found for the void that was embedded in a spatially 
flat Friedmann-Robertson-Walker (FRW) background when Hq was included, but they were 
removed by allowing the background to be hyperbolically curved. The authors in Ref. [16] 
found their models to fit to the full CMB spectrum and SNIa data, but the models suffered 
from a low Hq. Similar tensions between SNIa, Hq and CMB observations were encountered 
in Ref [5]. The LTB models can explain the kinematic Sunyaev-Zeldovich (kSZ) effect [17] 
only if SNIa, CMB and Hq fits are not considered [18-20]. The results from studies subjecting 
LTB models to observational data thus do not paint a clear picture of the pros and cons of 
the models. However, since different authors use different models, data sets and methods, 
the results are not necessarily contradictory. 

Our aim is to bring a different understanding to these results. We fix the luminosity dis¬ 
tance - redshift relation to be equivalent with that of the ACDM model. Then, we investigate 
which quantities differ most from their ACDM counterparts, and study observables which 
are dependent on the most diverging quantities. This type of LTB model has been studied 
before [21], but the issues at the apparent horizon (AH) have only recently been resolved [22] 
making it more appealing to study the model further.^ A disadvantage of this LTB model, 
where the observer is at the symmetry center and the bang time is homogeneous, is that it 
can mimic the distance - redshift relation only up to redshift z ~ 6.9 [22] after which further 
properties must be described by a model which no longer has the same distance - redshift 
relation as the ACDM model. 

LTB models are often analysed in terms of their specific relation to observations. This 
type of studies give us apprehension of the properties of the LTB models but not a full 
understanding about its failures and advantages. However, the fact that we have some sort 
of comprehension of the models capabilities, gives us a direction to go when generalising to 
more complex models; if we know where the LTB models appear to fail, it is wise to look 
carefully into these issues. 

In Section 2, we introduce the general properties of the LTB models relevant to this work 
and in Section 3, we derive the differential equations of the LTB model with homogeneous 
bang time, the ACDM distance - redshift relation and Hq in notation similar to that usually 

^The AH is the location where the derivative of the angular diameter distance with respect to the coordinate 
distance vanishes. 
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used in the ACDM model. In Section 4, we give an overview of how we solve the relevant 
differential equations and present the solutions; technical details of the procedure are given 
in Appendices A-D. In the following Section, 5, we compare quantities predicted by the 
obtained model with those predicted by the ACDM model. On the strength of these results, 
we continue by comparing these models with BAO (including Lyman-a forest) and CMB 
observations. We end the article with conclusions and a discussion in Section 6. 


2 The Lemaitre-Tolman-Bondi models 

The LTB models describe a spherically symmetric dust distribution. The metric is 


ds = —dt + 


1 + 2e(r)r2 

and the Einstein equations yield the evolution equation 

n, . 2M(r) . . 


dr^^ + R^{r, t){d9^ + sim Odcj)'^) 


( 2 . 1 ) 


( 2 . 2 ) 


where the functions M(r) and e(r) are time independent integration constants. The equation 
for the physical matter distribution is 


p{r,t) = 


2M'{r) 




(2.3) 


where k = SttG/c^. We denote partial derivatives with subscript r or t after commas. The 
solution of the former equation can be given in parametric form. 


R{r,t) = {cosh[77(r, t)] - 1} , 

. , , , .. , . [2e(r)r^]^/^ 

smh[r]ir,t)] - r]{r,t) = -- [t - tb{r)], 

where tb{r) is the bang time function. 

The redshift, z, and the coordinate distance, r, are related by 

= (1 + ~) 

dr y^l + 2e{r)r^ ’ 

where t should satisfy the null geodesic equation for incoming light rays 


(2.4) 

(2.5) 


dt 


R,r{r,t) 


dr -^1 + 2e{r)r‘^ 


( 2 . 6 ) 


(2.7) 


Using Eq. (2.4), the functions R^rir,t) and R^trir,t) can be given as functions of M{r), e(r), 
r and r](r,t). The dependence on the conformal time, r](r,t), can be eliminated using Eq. 
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(2.5), yielding 


R,r{r,t) = 


(2r^e'(r) + 4re(r)) ^ M'{r) 


2r‘^e(r) 


M{r) 


R{r, t) 


+ 

3 (2r^e'(r) + 4re(r)) {t — ti,{r)) 

M'{r){t - tb{r)) 

4r2e(r) 

M{r) 

1 

R^tr{r,t) = R,t{r,t) + 

2e[r) r 



- t'bir) 


Rtir,t), (2.8) 


M(r) 
' R{ry 


3 (2r^e'(r) + 4re(r)) (t(r) - 4(r)) M'{r){t{r) - th{r)) , 


4r^e(r) 


M(r) 


- t'bir) 


(2.9) 


The apparent horizon (AH) is dehned as the locus where dR[t{r),r]/dr = 0, where t(r) 
is given by the null geodesic equation (2.7) [23]. This definition reduces to the form 


R[r,t{r)] = 2M{r), (2-10) 

which can be seen by using the null geodesic equation (2.7). 

To ensure that the positive energy condition, p(r, t) > 0, is fulfilled, we demand M'{r) > 
0 and R^rit, r) > 0 throughout this paper (see Eq. (2.3)). Furthermore, to avoid shell crossing 
(SC) singularities, defined as surfaces where R^rir, t) = 0 and M'{r) / 0, inequalities t'^{r) < 0 
and 

^ [e(r)r2] > 0 (2.11) 

dr 

should always be satished [24]. 


3 The Lemaitre-Tolman-Bondi models mimicking the ACDM model 

We consider LTB models where the observer is located at the center, the bang time function 
is constant, and the Hubble constant and luminosity distance - redshift relation are equal 
to those in the ACDM model. Thus, a LTB model automatically gives hts identical to the 
corresponding ACDM model to Hq, SNIa and CMB dipole observations. 

Requiring that the LTB and the ACDM luminosity distances are equal is equivalent 
to requiring the angular diameter distances to be equal [21]. The ACDM angular diameter 
distance is given by 



(3.1) 


where is the “Friedmannian” Hubble parameter and the corresponding 

matter, curvature, and cosmological constant densities, respectively. In addition to a cos¬ 
mological constant, we have included a constant equation of state {p = wp) term The 
standard relation = 1 holds. 

The angular diameter distance for the LTB models is R[r(z),t{z)], thus imposing 

R[r{z),t{z)] = Da, 


- 4 - 


(3.2) 






















the LTB models have the same angular diameter distance as the corresponding ACDM mod¬ 
els. To ensure that the models behave as a homogeneous and isotropic models at early times, 
we set the bang time to a constant, which without loss of generality can be set to zero, i.e. 

tbir) = 0. (3.3) 


This choice immediately satisfies the condition < 0 required to avoid SC singularities, 
thus the only condition left to be satisfied to prevent SC singularities is Eq. (2.11). 

In order to write Eq. (2.2) in a more familiar form, we define 


M(r) 

e(r)r^ 

HR{r,t) 

HT{r,t) 

Ho{r) 


R,tr{r, t) 
R,rir,t) ’ 
R,t{r,t) 
R{r,t) ’ 

Hrir, to), 


(3.4) 


where Ro{r) = R(r,to) and to is the present time. Two Hubble functions are dehned above: 
Ht is the transverse Hubble function and Hr is the radial Hubble function. The geo¬ 
metrical mean of these Hubble functions in the three spatial directions is HLTB{r,t) = 

-I / Q 

[//ij(r, t)//j.(r, t)] ' . These definitions allow us to rewrite Eq. (2.2) as 


H^{r,t) 


Hiir) 


^M{r) 


( Mr) 

\R{r, t) 


3 

+ ^K{r) 


( Mr) 

\R{r, t) 


2 


(3.5) 


Because, at to for every r, HT{r,to) = Ho{r) and R{r,to) = Ro{r), the relation 


Hi^(r) = 1 - Quir) 


(3.6) 


must hold for every r at any given time. Using the gauge 

Mr) = r, 

Eq. (3.5) becomes 


H^{r,t) = Hlir) 


^uir) 


R{r, t) 


-|- (1 — Hm(?')) 


R{r,t) 


(3.7) 


(3.8) 


and at the homogeneous limit, where R{r,t) = a{t)r and Vt.M{r) and Ho{r) are constants, 
Eq. (3.8) reduces to the Eriedmann equation. 

Using the origin conditions given in Appendix C, we find HiTBiORo) = 77^(0, to) = 
HR{0,to) = Hq . Thus, to have appropriate conditions at the origin, the LTB models and 
the corresponding ACDM models must share the same local Hubble value, namely Hq . 

We will solve the system along null geodesics, where the derivative of R[r, t(r)] is 

d R 

~i~R = R r ~\~ R it r — R r — R t — ’ , ( 3 - 9 ) 

dr ’Vl + 2e(r)r2’ ^ ^ 
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and the second equality is obtained from Eq. (2.7). 

The system could be solved using Eqs. (2.6), (2.7), (3.2) and (3.9), which can be pre¬ 
sented as a set of differential equations for R, t, z and O. However, in order to improve the 
numerical stability of the solutions, we manipulate the set of differential equations into a 
more convenient form, see Appendix A. 

4 The solution of the system 


Age [GyrJ 




Figure 1. The LTB model with Hq = 67.1 and = 0.32. In the left panel, pieces of the functions, Hm 
(dashed orange), A — R/r (solid red) and « (dotted black), are depicted. Time, t, is plotted in the right 
panel. Fits appear perfect, and the redshift takes the value z « 6.9 at the end point r = rsc ~ 8.6. Colors 
available online. 

In Ref. [22], the system of equations in the special case where Hq and are fixed is 
studied. We study the system by allowing Hq and to vary. The set of equations solved in 
[22] differs from ours in two respects; the differential equations are constructed for different 
quantities and the gauges are chosen differently. We present here only an overview of the 
issues encountered solving the equations and direct the reader to Appendices A-D or [22] for 
more details. We encountered similar difficulties and numerical errors of the same order as 
those in [22]. 

The most challenging numerical task for numerically solving the system at hand is to 
find initial conditions so that the condition (2.10) holds at the AH. An algebraic relation 
between the initial and the AH conditions can not be found, hence one needs to choose 
initial conditions and numerically evaluate the system from the origin to the AH and check 
if the AH condition (2.10) holds. To preserve numerical accuracy at the origin and the AH, 
we use linear approximations near these points. 

The most accurate solution to the system is obtained by integrating the system from 
the vicinity of the origin to an intermediate radius (in Eigure 1 at r = 1.1 Gpc) and from 
the vicinity of the AH (in Eigure 1 at r = ~ 4.5 Gpc) to the same intermediate radius 

and then gluing these solutions together. Beyond the AH, we continued the integration of 
the system normally. The quantities t, R/r, z and H of the glued together solutions of the 
LTB model with Hq = 67.1 and = 0.32 are shown in Eigure 1. 

The plots in Eigure 1 ends at r = rsc ~ 8.6 Gpc, because the system violates the SG 
condition (2.11) there. Gonsequently, the LTB model can not describe a physically acceptable 
universe at this point when the AGDM luminosity distance with respect to redshift is imposed. 
However, the SG takes place at z ~ 6.9, where the AGDM model is dominated by matter. 
This allows us to continue the system further without adopting the luminosity distance of 
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the ACDM model. For our needs, it is sufficient to continue the system beyond the SC 
singularity by not breaking the SC condition and by demanding early time homogeneity. 
Thus, one should evaluate the physical matter density p{r, t) in the early times and make 
sure it is (approximately) independent of r. In this article, however, we use a weaker condition 
where II(r) is constant from the bang time to the drag epoch. This appears to be adequate 
for our needs. 

5 Comparison to ACDM model 

Our model is designed to inherit the luminosity distance from the ACDM model before the SC 
singularity. Hence, if some other quantities differ, observables dependent on these quantities 
may differ too. Therefore, first we study which quantities differ considerably between the 
models and then observables dependent on those quantities. For this, we use the LTB and 
the ACDM models with Hq = 67.1 and = 0.32. 

The difference between the age functions, t{z), is interesting, because the age of the 
universe has been used to constrain the LTB models in the literature. We find that the 
relative difference is approximately 6%, 12% and 17%, at redshifts 0, 1 and 6, respectively. 
This has been studied also in Refs. [16] and [25]. 

Instead of investigating quantities R^riz), R,t{z) and R^tr{z) separately, we investigate 
the transverse Hubble function and the radial Hubble function H^{z). The differences 

between H^{z), H^{z) and the ACDM Hubble function H^{z) are of similar size as the age 
function; at the origin, Hq = = 77^ and at redshift 6, the difference between 

H^{z) and H^{z) is approximately 16%, and between H^{z) and H'^{z) approximately 
18 %, see Figure 2. 

In the inhomogeneous models, the expansion depends on position, which is not the 
case in homogeneous and isotropic models. Hence, the differences between the expansion 
rates of the ACDM and the LTB models are expected. From the metric (2.1), one can 
see that the expansion of the universe in the radial and transverse directions is encoded in 
the functions R^rir,t) and R{r,t), respectively. An illustrative way to compare expansion 
rates is to evaluate how much R^r{r,t) and R{r,t) have changed from some early time te 
to the present. For this purpose, one should evaluate the ratios R^r[r{z),t{z)]/R^r[r{z),t(,], 
R[r{z),t{z)]/R[r{z),te\ and a[t{z)]/a{te) at different radial coordinates, where a{t) is the 
scale factor of the ACDM model. However, these ratios contain the inconvenient dependence 
on te, for which reason we compare the ratios 

^ R^r[riz2),te] 

RArizi),te] RAriz2),t{z2)]' 

T _ R[r{zi),t{zi)] R[r{z 2 ),te] , . 

R[T{zi),te] R[T{z 2 ),t{z 2 )y 
^ a[t(^i)] ajte) 
a{te) a[t{z 2 )]' 
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Figure 2. The LTB and the ACDM models with Hq = 67.1 and = 0.32. The left panel shows the 
radial Hubble function, H^{z) (dashed red), the transverse Hubble function, H^{z) (dotted blue), and the 
Hubble function of the ACDM model, H^{z) (solid blue). All three Hubble functions coincides at z = 0. 
The relative difference between H^{z) and H^{z) is approximately 6% and 16%, and the relative difference 
between (z) and H'^ (z) is approximately 13% and 18%, at redshifts 1 and 6, respectively. The right 
panel shows RLER’s, where TZ^, TZ^ and TZ^ correspond to dashed red, dotted blue and solid blue curves, 
respectively. The relative difference between TZ^ and TZ^ is approximatelly 43% and 47%, and the relative 
difference between TZ^ and 7Z^ is approximately 35 % and 45 %, at redshifts 1 and 6, respectively. Colors 
available online. 


We fix the second factor of the ratios by evaluating it at the origin, hence 




Rr[r{z),t{z)] 
Rr[riz),te] 
A[r{z),t{z)] 
A[r{z),te\ 

a[t{z)] a{te) _ 

a(te) a(0) 1 + z^ 


A{0,te), 

A{0,te), 

1 


(5.2) 

(5.3) 

(5.4) 


as R{r,t) = A{r,t)r, R^r{0,te) = A{0,te) and i?^r(0,to) = ^(0, to) = 1- The relative local 
expansion rates (RLER), TZ^, , and TZ^, represent the amount of local expansion from 
te to our null cone compared to the expansion at the origin. Going through the coordinate 
values from r(z = 0) to r(z = 6.9), we obtain information about how the universe has grown 
at different coordinate distances. RLER’s TZ^, TZ"^ and 7Z^ of the ACDM and the LTB models 
are drawn in Eigure 2. Because the amount of local expansion is scaled with the amount of 
expansion in the origin, TZ^, TZJ and TZ^ all take the unit value at r = 0. Moreover, TZ^, 
TZ"^ and TZ^ all decrease with radius, showing that the expansion of the LTB model has been 
fastest at the origin. There are two reasons for this; distant objects have had less time to 
expand and the surrounding matter distribution is denser at large r. 

While the dependence of TZ^ vanish completely, TZ^, TZ7 still explicitly include the 
parameter For practical calculations we set tg = 3 x 10“^ (corresponding to about 300 000 
years). For this value, the difference between p[r{z = 0),fe] and p[r{z = 6),te] is ~ 10“^ %. 


5.1 Baryon acoustic oscillations 

The results above indicate that the predictions for observables dependent on RLER’s and 
Hubble functions will differ considerably between the models. As we shall see, baryon acoustic 
oscillations are such observables. Here we compare the models using not only the low redshift 
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galaxy BAO measurements, but also baryon acoustic features obtained from the Lyman-a 
forest (LaF). 


5.1.1 Maximum likelihood analysis 

The characteristic size of baryon acoustic oscillations is given by 

dzBAo = > ( 5 - 5 ) 

V zbao J 

which with the line element (2.1) and the redshift equation (2.6) gives 
jLTB _ ( ^^phys, BApy ^phys,BAo(^ + ^BAo)R^'f'^ 

“ V RIao ZBAoRf^O 



at the lowest order of approximation. The subscript or superscript BAO indicates that the 
quantity is evaluated at zbao- 

To obtain BAO predictions for the LTB model, we need to evaluate the physical lengths 
at early times. As shown e.g. in Ref. [15], in the LTB model, the physical scales in transverse, 
Iphys radial, directions change as 


Iphysil^^ t) 

^physiol 


\ ^physiol ^e)i 


(5.7) 

(5.8) 


We require the LTB model to be almost homogeneous at early times. For this reason, we 
assume that the mean primordial perturbation have grown equally in all directions, i.e., 
te) ~ ^J/iys(^ 5 ^e)- Hence we can collectively denote physical lengths in each direction 
% Iphysi'i'jte)- Furthermore, it is also reasonable to expect the mean perturbation to have 
the same size everywhere, i.e., IphysiTBAO^te) ~ Iphysi'f’^'te) for any r?‘ Consequently, Eq. 
(5.6) can be written as 


d 


LTB 

Zbao 


( 1 + zbaoH^ao \^^^ ( 

Zbao ) \^R,r{'^BAO,te) 


1/3 


/ Rbao \ . 

\R{rBAO,R)) 


(5.9) 


where Eq. (3.2) and the definition for the radial Hubble function in Eq. (3.4) are employed. 
The ACDM counterpart of the above equation is 


d 


F 

^BAO 


( 1 + Zbao Rbao 


1/3 


1 + Ze 

1 + Zbao 


^phys (de)- 


(5.10) 


We present the results of the high redshift BAO (or LaF) surveys [30, 31] as 
and (1 + z)A9°^^ on Table 1, a manner convenient for this work. In spherically symmetric 
space-times 1/Az = /c and [(1 + 2;)A0°^'^]“^ = a_LT>^*‘^/r;^*'^, where superscript 

fid refers to a homogeneous fiducial model, to a comoving sound horizon at the drag epoch 

^The early times homogeneity assumption reduces here the degrees of freedom. For the inhomogeneous 
early time cosmology, however, it would be difficult to determine the functions and ifhysdy't)- On 

the other hand, we could use the BAO obsevations to determine these functions. 
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redshift 

dT 

AzObs 

(1 ± z)A0°'’* 

Ref. 

abbreviation 

0.106 

0.336 ±0.015 



[26] 


0.35 

0.113 ±0.0022 



[27] 


0.44 

0.0916±0.0071 



[28] 

low or 

0.57 

0.073 ±0.0012 



[29] 

low BAO 

0.60 

0.0726±0.0034 



[28] 


0.73 

0.0592±0.0032 



[28] 


2.34 


0.109 ±0.0033 


[30] 

LaFacR 

2.34 



0.08865±0.0051 

[30] 

LaFacT 

2.36 


0.111 ±0.0037 


[31] 

LaFccR 

2.36 



0.0926 ±0.0.003 

[31] 

LaFccT 


Table 1. The BAO observables used in this work. The observational values in Refs. [27, 29-31] are 
reported as l/do .35 = 8.88 ± 0.17, l/do .57 = 13.67 ± 0.22, l/A^ = 9.18 ± 0.28, [(1 + 2)A0]-i = 11.28 ± 0.65, 
I/A 2 = 9.00 ± 0.30 and [(1 + z)lS.6\~^ = 10.8 ± 0.4, respectively, but on the table we present their inverses. 
In the table, measurements are devided into low [26-29] and high [30, 31] redshift BAO measurements. 


and ay = and a± = A9'f^'^/A9. In homogeneous space-times also the relations 

1/Az = Hrd/c and [(1 + z)A9°'^^]~^ = hold, hence reporting the resnlts as Hrd/c and 

DA/f’d-, as is done in Refs. [30, 31], is merely a homogeneous interpretation of the results. 

We assume no correlation between the high and the low redshift observations, hence the 
combined for the low and a given LaF observation in the radial direction is 


{Xlow+LaF r) 


iXlow) 


2 

^LaFR 


and in the transverse direction is 


(5.11) 


Here 


(A/oiu+LoFt)^ = {Xlowf + 


[{1 + z)A9°’’^ - {1 + z)A9Y 


a 


LaFT 


Az^^^ 

Az^ 

A 9 LTB 

A9^ 


j^LcxF 

(1 “ 1 “ ^e^Hj^Qiplphysite), 

j^LaF 

(1 -|- ZLaF)Hp^prj^ ^ ^ s^lphys{te)i 

(1 “ 1 “ Ze^H]^aplp}iyg{teY 


(5.12) 


(5.13) 

(5.14) 

(5.15) 

(5.16) 


and the observed values are tabulated on Table 1. We study the following five cases separately: 

- low BAO 

- low BAO and LaF auto-correlation in the radial direction (LaFacR) 

- low BAO and LaF auto-correlation in the transverse direction (LaFacT) 

- low BAO and LaF cross-correlation in the radial direction (LyaFccR) 

- low BAO and LaF cross-correlation in the transverse direction (LyaFccT) 
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data sets 

model 

Xmm 

m,min 

a 

2a 

3cj 

1-p 


ACDM 

1.7 

0.23 

[0.11,0.39] 

[0.075,0.47] 

[0.048,0.54] 

0.21 


LTB 

2.6 

0.47 

[0.38,0.58] 

[0.36,0.64] 

[0.34,0.69] 

0.38 

low BAO 

ACDM 

1.7 

0.23 

[0.19,0.27] 

[0.18,0.29] 

[0.17,0.31] 

0.11 

+ LctFacR 

LTB 

15 

0.35 

[0.31,0.39] 

[0.30,0.41] 

[0.30,0.42] 

0.99 

low BAO 

ACDM 

1.7 

0.24 

[0.20,0.28] 

[0.18,0.30] 

[0.17,0.32] 

0.11 

+ LaFccR 

LTB 

13 

0.36 

[0.32,0.41] 

[0.31,0.42] 

[0.30,0.44] 

0.98 

low BAO 

ACDM 

2.1 

0.26 

[0.16,0.40] 

[0.13,0.48] 

[0.10,0.54] 

0.16 

+ LaFacT 

LTB 

2.9 

0.48 

[0.40,0.59] 

[0.38,0.65] 

[0.36,0.69] 

0.28 

low BAO 

ACDM 

4.0 

0.33 

[0.23,0.47] 

[0.19,0.53] 

[0.17,0.59] 

0.45 

+ LaFccT 

LTB 

4.8 

0.53 

[0.44,0.64] 

[0.42,0.69] 

[0.40,0.73] 

0.56 


Table 2. The minima of and corresponding for different data sets and models are tabulated. The 
minima appear to be independent on and the na-limits are given for OZi- The p-values are given as 1 — p 
on the last column. 


These five different cases are investigated separately because of the acknowledged discrepancy 
between LaF data and the ACDM model, which has raised suspicions of the validity of the 
measurements (see [33]). We expect these cases to reveal if discrepancy exists also in the 
void models. Moreover, we also expect to obtain some comprehension of the effects of the 
systematic errors. 

The analysis is executed so that for each and Hff we find the minimum Xmin 
both models by allowing Iphys to vary. We constrain parameters and lphys{te) so 

that ^ 

0.01 < < 0.85, 30 < < 80, 0 < lphys{te)- (5.17) 

We find the to be independent of for both models. Numerical analysis reveals that 
, and depend only on the combination {H^Y/'^lphysite), whereas df oc 

Az^ oc A6^ oc HQlpiiys{t(,) [see Eq. (5.10)]. Consequently, we can fix Hq using local Hubble 
observations solely, thus leaving us two free parameters to fit the BAO observations. The 
curves xfowi^m) both of the models are plotted in Figure 3, together with the latest 
SN constraints presented in [34].^ Both models can explain the used low BAO data with 
comparable minimum xfow Table 2) though the LTB model is in distinct conflict with 
the SN data. Also X^ow+LaFccRi^m) and xfow+LaFacA^m) curves are drawn in Figure 3. 
These specific xfow+LyaFi^m) curves were chosen, because they lay less stringent constraints 
to the models in radial and transverse directions compared with the alternatives (see Table 
2 ). 

For any of the five combinations of low and high BAO studied above, we can include 
the SN consraints [34] and calculate Xtot by 

xUO = x%Ao(i‘i) + p ig) 

’^obs 

^We were unable to execute the procedure described in Section 4 for < 0.2. Therefore, the was not 
calculated for the LTB model for these values. We identified the reason to be that <; q 2 requires 
fiA/fO) < 0. 

"^We use supernova constraints here for two reasons. We took the ACDM luminosity distance for the LTB 
model for the very reason the models fit equally well to the supernova observations. In addition, it is not 
evident how the Planck results [38] should be interpreted here. As they are obtained from the CMB sky, their 
feasibility constraining the inhomogeneous models that fits supernovae is questionable. 
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data sets 

model 

A,mzn 

m,min 

CJ 

2cj 

3(T 

l-p 

ABIC 

CAT 

ACDM 

2.2 

0.29 

[0.23,0.34] 

[0.21,0.36] 

[0.19,0.38] 

0.30 

^ n 

iUW “T olN 

LTB 

12 

0.37 

[0.33,0.40] 

[0.32,0.42] 

[0.31,0.44] 

0.98 


low -|- SN 

ACDM 

4.2 

0.25 

[0.22,0.29] 

[0.20,0.30] 

[0.19,0.31] 

0.48 

1Q 

-|- LaFacR 

LTB 

17 

0.34 

[0.31,0.37] 

[0.30,0.38] 

[0.29,0.39] 

1.00 

io 

low -|- SN 

ACDM 

3.5 

0.26 

[0.22,0.30] 

[0.21,0.31] 

[0.20,0.33] 

0.38 

1 9 

+ LaFccR 

LTB 

16 

0.34 

[0.32,0.37] 

[0.31,0.39] 

[0.30,0.40] 

0.99 


low -|- SN 

ACDM 

2.3 

0.29 

[0.23,0.34] 

[0.21,0.37] 

[0.20,0.39] 

0.19 

1 Q 

-|- LaFacT 

LTB 

15 

0.38 

[0.34, 0.42] 

[0.33,0.44] 

[0.32,0.45] 

0.99 

io 

low -|- SN 

ACDM 

4.2 

0.30 

[0.25,0.36] 

[0.23,0.38] 

[0.22,0.40] 

0.48 

1 Q 

-|- LaFccT 

LTB 

23 

0.40 

[0.36, 0.44] 

[0.35,0.46] 

[0.34,0.47] 

1.00 



Table 3. The minima of and corresponding for different data sets and models are tabulated. The 
minima appear to be independent on and the ncr-limits are given for The p-valne and ABIC are 
given on the two last colnmn, the p-value reported as 1 — p . The SN data is obtained from Ref. [34] . 


where we use = 0.295 ± 0.034 obtained from [34]. The results are presented on Table 3. 

We use Bayesian information criteria (BIC) and p-value to asses the strength of the 
models. Because both models have equal amount free parameters and the used errors are 
Gaussian, the difference in the BIC values of the models is 

ABIC = (5.19) 

A difference in BIC of 2 is considered positive evidence against the model with the higher 
BIC, in this case higher x^ corresponds to higher BIC, while a difference in BIC of 6 (or 
more) is considered strong evidence [35]. For further information, we refer to [36] and [37]. 

All the evidence obtained here using the maximum likelihood analysis favours the ACDM 
model. Especially for the combined SN and BAO data, the ABIC values exhibit strong evi¬ 
dence against the LTB model and p-values rule out the LTB model at least at 98% confidence 
(see Table 3). Figure 3 illustrates the effect of LaF data for both models. The discrepancy 
between SN, low BAO and high BAO is present for both models. For the LTB model, the 
transversal LaF data appears to be consistent with low BAO, but in severe conflict with SN 
data, whereas the radial LaF data appears not to conflict much with SN data, but is in a clear 
contradiction with low BAO observations. The weaker constraining curves drawn in Figure 
3 demonstrate the situation well on a qualitative level and using the stringent constraining 
curves would not modify these conclusions. Therefore, we conclude the systematic errors in 
LaF observations do not play a significant role here. 


5.1.2 Scale independent analysis 

In addition to the x^ analysis, we compare the models using a method independent on Iphys- 
For this analysis, we choose to use LaF data for which the values of the models are 
comparable and value acceptable in the light of SN data. For these reasons, we choose 
to use low BAO and LaFccR data and fix 0^ = 0.32. Taking the ratio dzjAz eliminates the 
dependence on Iphys and at redshifts zi and Z 2 we find 


dLTB 

Az^'^^ 




[(1 -h Z2)H2] 


nf) 


1/3 



2/3 


(5.20) 
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Figure 3. The left panel: Different x^(^m) curves for different models and data sets are presented. The 
black vertical line in the middle of the vertical contours indicates the preferred according to the SN data 
[34] , whereas the vertical contours indicates 1-3 a deviations from the best fit value. See Table 2 for numerical 
values. The right panel: The (black) points and bars represent observations and their Icr errors, the red 
solid curve represents the ACDM model’s prediction and the dashed blue curve represent the LTB model’s 
prediction. For both models, = 0.32 and Az is chosen to correspond LaFccR. Colors available online. 


The ACDM counterpart of the above equation is 



(5.21) 


The difference between Eqs. (5.20) and (5.21) is caused by the Hubble functions and the 
relative local expansion rates of the different models. The ratios (i^^/Az 2 , where Az 2 = 
AzLaFccR, are plotted in Figure 3. The models appear to fit equally well to the observations, 
but calculating the variance, 



which for LTB is ~ 0.11 and for ACDM is ~ 0.19, reveals that the LTB model gives slightly a 
better fit. This is in minor disagreement with the maximum likelihood analysis, which gives 
a slightly lower value for the ACDM than for the LTB, as can be seen from Figure 3. 
The small disagreement can be understood by considering the nature of the ratio analysis. 
Because the variance does not take standard deviations into account, is fixed and Iphys 
is eliminated (unlike in the maximum likelihood analysis where they were optimised), is 
variance (5.22) merely a measure how much data points differ from fixed models. From 
Table 2 can be seen that the difference between of low BAO + LaFccR and the fixed 

Dm = 0.32 is 0.4 for the LTB and 0.8 for the ACDM model. Even though the values in 
the table correspond to maximum likelihood analysis, we expect the corresponding 
values for the variance to be in the vicinity, thus explaining the small difference in the results 
of two analysis methods. 

The ratio analysis demonstrates how the differences between the Hubble functions and 
RLER’s of the two models are inherited by the measurable BAO quantities. Figure 3 demon¬ 
strates how the difference in the predictions of the BAO observables of the models increases 
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with respect to redshift, as expected (see Figure 2).® At the same time, the errors appear to 
grow with redshift, hence the constraining power of the wide range of the observations is not 
strong with this method. 

We deduce the results obtained here are mainly caused by the RLER’s for two reasons. 
The RLER’s differ more than the Hubble functions at the redshift range of the observations. 
Moreover, in Eqs. (5.20) and (5.21) is effectively proportional to the power of —2/3 

of the Hubble functions, whereas they are effectively linearly dependent on RLER’s. 

Both of the two analyses showed that BAO observables do inherit the growing con¬ 
straining power with respect to redshift from the Hubble functions and LRER’s. Minimum 
likelihood analysis clearly favours the ACDM model, but demonstrated the models to be 
comparable if LaE data in the radial direction is employed and is close to the value 
favoured by Planck [33]. Analysing this specific case with a scale independent method, we 
found the LTB model to be slightly favoured. However, it is not evident how homogeneously 
interpreted Plank results should be received here. Moreover, we expect that the scale inde¬ 
pendent analysis using transverse LaE data would yield a contrary result. On the strength 
of these results, we conclude that the ACDM model fits combined BAO and SN observations 
better than the LTB model. 


5.2 Cosmic microwave background 

The angular diameter of the decoupling surface d* is in the LTB model given by 

II 


e* = 


^physir*X) 


R{r*,t* 


(5.23) 


In practical calculations, we take the decoupling time t* to correspond to a redshift z = 1089. 
The latest Planck results [38] give 9* = (1.04147 ±0.00062) x 10“^. Since at early times, the 
mean physical lengths are assumed to be the same at every r, Eq. (5.9) can also be written 


1/3 


d 


■LTB ^ / 1 ± ZB AO \ 

ZBAO V^,r-(0,td) 


'T^R 

'^BAO 


1/3 


/ 


BAO 




2/3 


^physi^dRd) ^ (^’2^) 


where the index d refers to drag epoch, for which we use the redshift z = 1020. At the 
homogeneous era 


d +d\ rr'^+Ts 


lphys{r^,t'^) 


dr 




dr 1 z* r^ 


Iphysir* ,t*) R^r{r*,t*) dr a(t*) 


r*+r* 


dr l + z^ r* 


1090 

1021 


(1.0188 ±0.0068), 
(5.25) 

where the numerical value for r//r* is taken from Ref. [33]. Combining (5.23)-(5.25), we can 
write 


dLTB 

^zbao 


/ 1 ± ZbAO ^bAO \ I ''^BAO \ / ''^BAO \ I + Z 


9* V {DBAOy 


) ( 


K ) 


V 


1 + z'^ r* 


R(r*,C), (5.26) 


since 


R,r{0,td) = ^(OTd) = = T^d- ACDM counterpart equation is 


^ZBAO ^ f 1 ± ^BAO \ 

9* V zbao J 


1/3 


'^BAoi'i^ + ^*HD*A- 


(5.27) 


®Note that the curves of the models do not converge when z grows larger than presented in Figure 

3. Instead, they cross at 2 ~ 0.85 and d^^/Aa is 0.27 for the ACDM and 0.26 for the LTB at 2 = 2.34. 
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Eqs. (5.26) and (5.27) include quantities dependent on the BAO, the drag and the 
decoupling epochs. By choosing how to continue the void profile function beyond the SC 
singularity (the tail) one fixes the values of the quantities dependent on the drag and 
the decoupling epochs. On the other hand, we have fixed = AziaFccR- Therefore, 
oc and dfj/Azf’ oc dzgAo/^*' (^-26) and (5.27) predict the BAO 

observables as rescalings of Eqs. (5.20) and (5.21). This is illustrated in Figures 3 and 4. 
This implies that regardless how the void profile is continued beyond the SC, the LTB model 
can hardly describe the wide range of BAO observables better than in the case where only 
BAO was investigated. 

Despite of the rescaling effect noted above, it is not evident that optimal tails are natural 
for the LTB model. Therefore, we consider several tails (illustrated in Figure 4) to be able to 
compare the models. The tails are chosen requiring that the SC condition is not violated and 
that the system is approximately homogeneous at z > 1000. The latter condition is satisfied 
by constant tails for z > 1020. However, this condition does not guarantee the universe to be 
homogeneous at z ~ 1020, but merely enables it (see Section 4). The tails in the left panel 
of Figure 4 are constructed as follows: The blue tail is rising rapidly and changing to a 
constant before the SC singularity is violated. The orange curve D 2 is tailored to rise slightly 
slower than the blue curve and changes to a constant before the SC singularity occurs. The 
red tail D 3 is rising to encounter the SC singularity and the homogeneity limit almost at 
the same time. The brown tail D 4 rises slowly before it reaches the homogeneity limit. The 
green tail D 5 is set to be a constant at first, ascent after a while and turn to a constant again 
before the SC condition is violated. The purple curve De is set to a constant from z ss 6.9. 

In the right panel of Figure 4 we have plotted the LTB and the ACDM predictions for 
the ratio dzjO* using different void profiles together with the observed values. The models 
in question are characterised by Hq = 67.1 and = 0.32. As can be seen from the figure, 
changing the void tail does not seem to have great impact on the predicted ratios. Even 
though the LTB models ht well inside the error bars, they are characterised by the feature 
that they do not fit observations as well as the ACDM model. The LTB models can be 
arranged from the best fit to the worst fit. It appears that higher Qm maximum value and 
a steeper incline gives a better fit to the data. The problem of the linear piecewise functions 
is that they can not be as steep and rise as high as needed without violating either the shell 
crossing condition ( 2 . 11 ) or the homogeneity assumption. 

We also applied more complex functions to describe the tail and considered the possi¬ 
bility of it starting at smaller redshifts. Evidently, this causes the luminosity distance of the 
LTB model to begin to deviate from the ACDM model at lower redshifts. Although tails with 
a steep ascent to high maximum values made the fit better, they did not differ considerably 
from the piecewise tails. 

The difference between Eqs. (5.26) and (5.27) is caused by the Hubble functions, the 
RLER’s and the angular diameter distances after the SC. Although both of the equations 
include the factor r^/r*, a somewhat reasonable assumption is that it does not take the same 
value in both models. After all, r^/r* ~ 1 implies assuming the baryon fraction is the same 
at the BAO, drag and decoupling epochs. With a variable baryon fraction one would obtain 

^ = c(1.0188± 0.0068), (5.28) 

where c is a constant, determined by the amounts of baryonic and dark matter at and 
fBAO- can choose this constant to give the best possible ht for the LTB model, which 
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Figure 4. The LTB and the ACDM models with Hq = 67.1 and = 0.32. In the left panel several tails 
are plotted. In the right panel, (black) points and bars are observed d^jO* values and their Icr error bars, the 
solid blue curve represents the ACDM prediction and non-solid colored curves represent the LTB models with 
tails, the curve color and style indicates the correspondning left panel tails. The LTB curves can be arranged 
from the worst to the best fit as follows: fie, fls, fli, 1^4, 0.2, and fls. Colors available online. 

corresponds to c Ri 0.7. Nevertheless, as we showed earlier, and are 

just rescaled and d^^^^/6*, respectively, thus the combined BAO and CMB set is 

at the best at the same level as the ht obtained using BAO observations only. Therefore, we 
conclude that the inclusion of CMB observations emphasizes the superiority of the ACDM 
model over the LTB model. 

Note that the LTB models are dust only solutions of the Einstein equations. Conse¬ 
quently, the effects of radiation are neglected here and for the reason the treatment of the 
CMB is incomplete. However, our analysis starts ai z = 1089, when the radiation density 
was still subdominant while not completely negligible. Therefore, our studies of the CMB 
are pointing how difficult it is for the LTB void to mimic different features of the ACDM 
model at wide redshift range. Moreover, because the ACDM model is in good concordance 
with the hrst CMB peak, we are conhdent our conclusion of the inclusion of the CMB data 
holds even if the effects of pressure would be considered. For discussion about the dynamical 
effects of radiation in inhomogeneous models, see [39]. 

6 Conclusions and discussion 

We have compared the ACDM and the LTB models which have the same distance - redshift 
relation and local Hubble parameter value. In addition, the LTB model is forced to have 
a homogeneous bang time and to be homogeneous between the bang time and the drag 
epoch. The relative local expansion rates appeared to have the largest difference, but the 
Hubble functions showed considerable deviation too. The difference between the ACDM and 
the LTB models is larger the wider the redshift range covered. On the strength of this, we 
investigated baryon acoustic oscillation and cosmic microwave background observations. The 
results favour the ACDM model. Furthermore, the LTB models exhibits conflict between SN 
and BAO observations (with or without LaF data). They suggest, however, that to obtain 
comparable fits for SNIa, Hq, CMB, and BAO, one needs to be able to control the LTB 
model’s luminosity distance, Hubble functions, and relative local expansion rates. To achieve 
this, the model needs more degrees of freedom than those employed in this paper. 
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The analysis of low redshift BAO observations slightly favoured the ACDM to the LTB 
model. Using the maximum likelihood analysis by allowing and Hq to vary, we noticed 
that the minimum value is independent on Hq for both models. For low BAO observations 
Xmm ACDM model is 1.7 and for the LTB model is 2.6. Even though these values are 

close to each other, for the LTB model it was obtained when the mimicked ACDM model has 
= 0.47, which is in conflict with the latest SN [34] and Planck [38] data. Thus, although 
both models can comparably explain low BAO observations, the combined low BAO and SN 
data clearly favours the ACDM model. This deduction is endorsed by p-value, which states 
the LTB model to be ruled out by 98% confidence and Bayesian information criteria stating 
there is strong evidence against the LTB model. 

We also investigated the effects of the baryonic features of the Lyman-a forest. For 
the combined set of low redshift BAO and LaF data in the radial direction, the LTB model 
predicts the set distinctly worse than low BAO only, but the tension to SN does not grow 
appreciably. The effect on the ACDM model is reverse as adding the radial LaF observations 
have almost negligible effect on explaining the BAO, but tension to SN grows drastically. 
For the combined BAO and supernova data p and BIC values clearly favour the ACDM 
model, even though the results are not complementing it either. The weakest constraint 
from p-values rule the ACDM out by 38% confidence and the LTB out by 98% confidence. 
Moreover, the BIC value indicates strong evidence against the LTB model. Furthermore, if 
low BAO data is combined with the transverse LaF data, the results are different from the 
radial case. For the LTB model, the transversal LaF data appears to be consistent with low 
BAO, but in severe conflict with SN data, whereas these three data sets flatter the ACDM 
model. The weakest constraint from p-values rule the ACDM out by 19% confidence and 
the LTB out by 99% confidence and the BIC value indicates strong evidence against the 
LTB model. In conclusion, the ACDM model fits BAO (with or without LaF data) and SN 
observations better than the LTB model. 

In addition, we compared how LTB and ACDM compare to BAO and first CMB peak 
observations. The LTB model can not mimic the ACDM luminosity distance - redshift 
relation at the shell crossing singularity, of which location depends on the parameters of the 
mimicked ACDM model; for parameters Hq = 67.1 and = 0.32 the location is at z ~ 6.9. 
For this reason, we modified the void profile beyond z ~ 6.9 for this particular model, and 
compared it with observations. We studied different tails, but no significant improvement was 
found. Only if the mean sound horizon depends on the position, some remedy was achieved. 
This possibility is related to the assumption of inhomogeneous early time, indicating our 
assumption of its homogeneity is too restrictive. Another deficiency in our approach is that 
the LTB model does not include pressure, thus leaving our CMB analysis incomplete. For 
these reasons, we conclude that under the assumptions made in this study, combining CMB 
to the BAO observations emphasises the better fit of the ACDM model. However, due to the 
concordance of the ACDM model and the CMB observations, we expect this to hold quite 
generally. 

Explaining the local Hubble value observations does not cause difficulties to the LTB 
model studied here, unlike in some other studies. The reason is our parametrisation, which 
allows us to satisfy the Hq observations independently on other observations. This is not 
the case e.g. in [14], where the other observables fix all free parameters and Hq is then to 
be calculated. In [15] the parameterisation allows freely to choose Hq for the model, but 
the authors find observably viable values to correspond to too young universe. This suggests 
that the issues confronted in explaining Hq observations could vanish for some LTB models 
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by reparametrisation. 

Our results agree with studies [5, 15, 16], but appears to be in conflict with [13] and [14], 
The LTB void models were found to have comparable fit to the ACDM model using SNIa, 
first CMB peak and BAO observations in [13] and in addition to Hq observations in [14], 
even though in the latter case, difficulties creating high enough Hq value were encountered. 
We identify the factors explaining these differences. In these articles, the sound horizon 
was implicitly allowed to vary with distances. Moreover, in both papers, the redshift range 
of the BAO observations was narrower even than the low BAO data in this article. We 
demonstrated how the deviation of the BAO features between the two models grows with 
respect to redshift. Because the trusted low BAO observations are in concordance with the 
ACDM model, it appears to be justified to say that the wider redshift range imply worse 
concordance with the LTB model and low BAO observations. This is backed by Ref. [15], 
where as wide redshift range of BAO observations was utilized as with low BAO here and 
the fit between observations and the models was far worse. 

We believe that we have found the reasons why LTB models constantly fails to explain 
Hq, SNIa, BAO and CMB observations simultaneously. Combining SNIa observations, that 
determine the luminosity distance - redshift relation, with BAO observations, which are heav¬ 
ily dependent on the local expansion rates, tension is manifested between the model and the 
observations. Extending the range from where the observations dependent on the local ex¬ 
pansion rate and the Hubble function are obtained, such as including CMB and local Hubble 
parameter in the range, the tension grows. This suggests that to obtain observationally more 
viable models, one needs to include more degrees of freedom by renouncing the homogeneous 
early universe hypothesis, the homogeneity of the bang time, the negligible effects of pressure 
or the spherical symmetry assumption. 
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A Enhancing the accuracy of the solutions 

To increase the numerical accuracy of the solutions of the system, we modify the system of 
differential equations (2.6), (2.7), (3.2) and (3.9). 
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Let us begin by eliminating the integral from Eq. (3.1) first as 


arcsinh I^DaH^{I + — 


dz 


+ Z)3 + Of (1 + Z)2 + Of (1 + Z)^ + Of 

(A.l) 

and then differentiate both sides with respect r, yielding 


aoR'{r) + ai2:'(r) = 0, 


(A.2) 


where oq and ai are auxiliary and A(r) = R{r)/r. We define a set of auxiliary functions Oj, 
z = 0, ...9, given below. 


ao — 


ai = 


Hq {z{r) + 1) 


)20fr2A(r)2(z(r) + 1)2 + 1 ' 
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(A.IO) 

(A.ll) 

(A.12) 


We represent Eqs. (2.6), (2.7) and (3.9) using auxiliary functions Oj. The common factors 
for all Oj are that they all take a non-zero value at the origin and they all depend only on 
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^Mir), z(r), t{r), A{r) and r. H^ir) is actually a function of r through only, 
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H^{r) = F, 


00 


1 — n(r) 


0(r)arcsinh 

(1 -Q(r))3/2 


(A. 13 ) 


With the help of the definitions ( 3 . 4 ) and the gauge equation ( 3 . 7 ), Eqs. ( 2 . 6 ), ( 2 . 7 ) and ( 3 . 9 ) 
now becomes 


_ «6 + {ajn'j^ir) + 08t'(,(r))r 
z — , 

05 

^ 02 + (a3^M(r) + a4t^r))r 
R'{r) = 09(02 + (o3n'(r) + 044(r))r]. 


(A. 14 ) 

(A. 15 ) 

(A. 16 ) 


Taking the homogeneous bang time to be ti, = 0 , using the dehnition R{r) 
rearranging (A.2), (A. 14 ), (A. 15 ) and (A. 16 ) we obtain 

= A{r)r and 

// N ao(“ 3«6 — 0207)09 

^ in - - - ^ — - -- , 

0107 -|- 09030509 

(A. 17 ) 

,// N 010309 — OIO2O7 

Hn = ^ 2 ’ 

030507 -|- 00030^09 

(A. 18 ) 

, 0i( — 0306 + 0207)09 

A (r) = 02 + 

QjICLy + CLQCL^CL^CLg 

(A. 19 ) 

! OlOg + 09020509 

— 

0i07r + CLQCL^a^CLgr 

(A.20) 


This set of differential equations gives considerably more accurate solution than equations 
( 2 . 6 ), ( 2 . 7 ), ( 3 . 2 ) and ( 3 . 9 ). 


B Solving the differential equations 

We follow the procedure to solve the system described in [ 22 ], but in a different gauge. We 
choose a gauge which at the homogeneous limit approaches the conventional ACDM model, 
making comparisons more straightforward. Because the differential equations (A. 17 )-(A. 20 ) 
are not well-behaved close to the origin and the AH, we solve the system in pieces. In the 
vicinities of the origin and the AH, the system will be solved using linear approximations. 

In this appendix we reconstruct the luminosity distance of the ACDM model given in 
[22], where the parameters of the homogeneous universe model are: 

ni = 0 . 32 , = 0, = 0, = 67 . 1 . (B.l) 

The results are obtained similarly for different parameters. 

B.l Solving the system starting from the origin towards the apparent horizon 

Because the set (A. 17 )-(A. 20 ) is singular at the origin, we use the linearized equations of 
the system from the origin to r = 10 “®. From there on we use Eqs. (A. 17 )-(A. 20 ). The 


- 22 














Figure 5 . The numerical evaluations for A, z, and t of the differential equations (A.17)-(A.20) are 
plotted. In the left panel, the evaluation begins at r = 0 and ends at r « 1.56545 x 10“^, where the 
numerical procedure hits a singularity. In the right panel, the evaluation begins at r = rAH and ends at 
r ~ 6.84513 x 10“®, where the numerical procedure hits a singularity. Colors available online. 


initial conditions are solved in detail in Appendix C, which with our parametrisation yield 
the numerical values (with six digits accuracy) 

Om(0) = 0.109621, t(0) = 0.0132964, A(0) = 1, z(0) = 0, , , 

Q'^(O) = 28.5060, t'(0) = -1, A'(0) = -67.1000, z'(0) = ^ ^ 

Using these initial conditions, the solution of the differential equations (A.17)-(A.20) and 
its linearization are plotted in Figure 5. The calculation stops in the vicinity of the AH, at 
r ~ 0.016, where the numerical methods report a singularity. 

B.2 Solving the system starting from the apparent horizon towards the origin 

We use the linearized equations to calculate the quantities from to — 10“®, where 
vah is the comoving coordinate at the AH, and use full Eqs. (A.17)-(A.20) for smaller r. 
The relevant initial conditions are presented in Appendix D and their numerical values with 
six digits accuracy are 

n M, AH = 0.518126, tAH = 0.00362764:, = 0.389397, z ah = 1-58243, 

^'m,ah = 23.4305, = -0.360968, = -25.3158, z'^^ = 170.758, ^ 

where the AH lies at r = tah = 0.0153816 with six digit precision. The solution is depicted 
in Figure 5, where functions ^(^)j ^(^)) z{r) are plotted. The numerical method 

hits a singularity at the origin. 

The solutions from the origin to the AH and from the AH to the origin are glued together 
at r = 0.0011, where the absolute difference of each quantity between both solutions is less 
than 2 x 10“® and the relative differences are all less than 0.07%. 

B.3 Solving the system beyond the apparent horizon 

To go beyond the AH, we solve Eqs. (A.17)-(A.20) starting from r = tah + 10“® towards 
infinity, nsing (B.3) as the initial conditions. The system hits the SC singularity at r = 
Vsc ~ 0.0285487, corresponding to a redshift of z = 6.93850, which is the same value reported 
in [22] within a precision of three digits. In Figure 6, the solutions is drawn in the range 
{rAH: 'i^sc) and its linear approximation is drawn at the vicinity of the AH. 


- 23 










7 


z(r),, • 


n(r) 


A(r) 


0.018 0.020 0.022 0.024 0.026 0.028 


Figure 6. The numerical solutions for Hm, A, z, and t (dashed blue, overlapping with the horizontal axis) 
of DE’s (A.17)-(A.20) are plotted. The evaluation begins at r = tah and ends at r = Vac, where the system 
hits shell crossing singularity. Colors available online. 


C The origin conditions 

To obtain a physically acceptable solution of the set (A.17)-(A.20), it has to satisfy the 
relation (2.10). It is straightforward to see that the system is regular at the origin, if 
liuir^Q R{r,t) = Ro{r) = r, limr^oM{r) oc r^, We can utilise the analysis of [22], where 
it was found that at the origin 


X^ + kX - 2 MqH^ = 0, (C.l) 

where X = Ao{0)Hq, k = —2e(0) and Mq = M(r)/r^|r=o- Using X, k, and Mq, from Eqs. 
(2.4) and (2.5) we find the present age of the universe, 


m = 


Mo 


(_fc)3/2 


-k X X , {-k X ^ 

+ 1 — 1 — arccosh —-^ + 1 




MoX 


from Eqs. (C.l) and (3.4) we get X^[Mo(0)^ — = 0, and Eqs. (3.4) reveal 

^MiO) = 

It is straightforward to see that the values 


(C.2) 


(C.3) 


X = 67.1, Mo = 100, X = 49.6539 and k = -2195.24. 


(C.4) 


satisfy Eq. (C.l) and the system solved from the initial conditions obtained from these values 
agrees with the condition (2.10) at high precision. 

Because we have defined R{r,t) = A{r,t)r, the gauge (3.7) dictates Ao(0) = 1. Using 
the L’Hopital rule, one finds 


lim 

r^O 


R{r) 


y /o^ [^m(l + ^)^ + + ^)^ + + ^a] ^ dz 

r“o Xr 
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r—^0 


dz/dr 
ffF 
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(C.5) 
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thus 


dz 

dr 


= H, 


r=0 


0 • 


(C.6) 


To find the value of dt/dr in the origin, we combine Eqs. (2.7) and (3.9) and use the chain 
rule, when 


R,r{r,t) = 1 - 


t) 


-1 


d[{r{z),t{z)] dz 


+ 2e(r)r2 


dz 


dr 


(C.7) 


On the other hand, differentiating both sides of (3.1) with respect to z, yields 


dR[r{z),t{z)] 

dz 


—z 


(1 + z) Hi (1 + z)2j (1 + ^)3 + 

1 (—311^ — 4)z 


+ 0(z^) 


■”o 


+ 


2Hi 


+ Oiz^ 


and, because 77^(0, to) = Hq , Eq- (3-8) yields 

R,tir,t) = Hir + 0{z^). 
Inserting Eqs. (C.6)-(C.9) into Eq. (2.7) gives 

= - 1 . 


dt 

dr 


r=0 


(C.8) 


(C.9) 


(C.IO) 


The values of dVt{r)/dr and dA{r)/dr at the origin we obtained by Taylor expansion of Eqs. 
(A.17)-(A.20) and finding the limit r —)■ 0. The result is 


d^uir) _ _ {Hjf (flM(O) - 1) f2M(0) (30^ - - 2) _ 

- 1) ^m(O) + Hi (flM(O) + 2)) {Hito (Om(0) + 2) - 2) ’ 

(C.ll) 

dA{r) _ f Hj (312^ + OM(r) + 2) (g^tp (0^(0 - 1) ^M{r) - l) 

I {^SV) - 1) OM(r) + Hi (OM(r) + 2)) to (f^M(r) + 2) - 2) 
~^dnM{r) ~ ^M{f) + + 2 ) {3Qrn + + 2 ) 

2 (2g^ (f2M(r) - 1) ^M{r) + Hi {^Mir) + 2)) {Hito (OM(r) + 2) - 2)' 

(C.12) 


Finally, we note that the Eqs. (2.6), (2.7), (C.6) and (C.IO) yield H^{0,to) = Hi, and, 
as 77^(0, to) = Hi, the mean Hubble value of the LTB model in the origin is exactly same 
as in the ACDM model. 


D Apparent horizon conditions 

In this appendix we mimic the luminosity distance of the ACDM model given in [22] (the 
relevant parameters are given in (B.l)). Similarly the results are obtained for different 
parameters. 
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The redshift and the angular diameter distance at the AH are found where the derivative 
of the angular diameter distance (3.1) with respect to 2 vanishes. The corresponding values 
are z = zah ~ 1.58243 and R[t{rAH)-,'fAH] ~ 0.00598955. Numerical methods fail to evaluate 
Eqs. (A.17)-(A.20) properly at the AH; this can be traced to Eq. (A.16). By expanding og in 
the vicinity of the AH, taking into account the AH conditions dR\r,t{r)\/dr = 0 and (2.10), 
it can be written as ag = [012 + {r)]{r — tah) + O [(r — rAn)'^], where 012 and 013 are 

given in (D.4) and (D.5). On the other hand, we can solve R'{r) from Eq. (A.2) and expand 
it around the AH, yielding R'{r) = aioaiiz'(r)^(r — tah) + O [{r — vah)"^], where aig and 
ail are given in (D.2) and (D.3). Equating the expressions for R'{r) we obtain 

aioauz'{rf = [012 + ai3n'(r)][o2 + (a30'(r) + ait'f^{r))r]. (D.l) 


The values for A'{r), f2^(r), z'(r) and t'{r) at the AH can now be solved from the equations 
(2.10), (A.14), (A.15) and (D.l) for the linear expansion around the AH. Eq. (D.l) shows 
that there are two solutions; we choose the one where t'{r) < 0 . 

The functions aio, ••■,013 are given by 
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ai 2 
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(D.2) 

nF + if^ 

(D.3) 

(D.4) 

(D.5) 


E The units 

In this article we use units where the speed of light c = 1 and cH^ = 67.1 km/(s x Mpc). 
This yields, within three digit precision, 

9.78 X lO^Gy = 1, (E.l) 

3.00 X 10^ Gpc = 1. (E.2) 

Thus, the age of the LTB model is 13.0 Gy (when the Friedmannian parameters (B.l) are 
assumed) and the cases where homogeneity is imposed at the drag epoch, void sizes are ~ 10 
Gpc. 
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